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ABSTRACT 

We revisit the pro ximity effect produ ced by QSOs at redshifts 2.1 — 3.3 applying 



t-^ . We revisit trie pro ximity ettect produ ced by QbUs at redsnilts z.i — 6.6 applying 

the FLO approach (jSaitta et al J l2007h to a sample of ~ 6300 Lya lines fitted in 21 
high resolution, high signal-to-noise spectra. This new technique allows to recover the 
hydrogen density field from the Hi column densities of the lines in the Lya forest, 
on the basis of simple assumptions on the physical state of the gas. To minimize the 
systematic uncertainties that could affect the density recovering in the QSO vicinity, 
we carefully determined the redshifts of the QSOs in our sample and modelled in 
detail their spectra to compute the corresponding ionising fluxes. The mean density 
. field obtained from the observed spectra shows a significant over-density in the region 

within 4 proper Mpc from the QSO position, confirming that QSOs are hosted in high 
density peaks. The absolute value of p/ (p) for the peak is uncertain by a factor of ~ 3, 
depending on the assumed QSO spectral slope and the minimum H I column density 
detectable in the spectra. We do not confirm the presence of a significant over-density 
extending to separations of ~ 15 proper Mpc from the QSO, claimed in previous works 
at redshifts (z) ~ 2.5 and 3.8. Our best guess for the UV background ionisation rate 
based on the IGM mean density recovered by FLO is Ttjvb — 10 -12 s _1 . However, 
values of Ttjvb — 3 x 10~ 12 s" 1 could be viable if an inverted temperature-density 
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' relation with index a ~ —0.5 is adopted. 



Key words: intergalactic medium, quasars: absorption lines, cosmology: observa- 
tions, large-scale structure of Universe 



1 INTRODUCTION 

The ultraviolet radiation emitted by quasars (QSOs) is con- 
sidered the dominant source of ionisation of the intergalactic 
medium (IGM) at redshifts 2-4. Most of the absorption lines 
seen blue-ward of the Lya emission in QSO spectra (the so- 
called Lya forest) are ascribed to fluctuations in the low to 
intermediate density IGM (see iMeiksin 2007, for a recent 
review). As a consequence, Lya lines can be used as probes 



* Based on observations collected at the European Southern Ob- 
servatory Very Large Telescope, Cerro Paranal, Chile — Programs 
166.A-0106(A) and during commissioning and science verification 
of UVES 



of the properties and redshift evolution of the UV ionising 
background. 



Observations show that the number density of Lya lines 
increases with redshift, but within single QSO spectra the 
number density of Lya lines decreases as the redshift ap- 
proaches ithe_QSO_eir^sionredshift. This effect was first no- 
tice d by Carswell et alj ( 19821') and co n firmed by late r stud - 
ies (|Murdoch et alj|l98d : lTvtlerlll987l ). [Baitlik et all (Il988l ) 
called this deficiency of Lya absorptions near the back- 
ground QSO "proximity effect" and attributed it to the in- 
creased ionisation of the Lya clouds near the QSO due to its 
ionising flux. They used the proximity effect in 19 low resolu- 
tion QSO spectra to estimate the intensity of the ultraviolet 
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background radiation (UVB) at the Lyman limilQ, Jll, for 
which they found the value log Jll — —21.0 ±0.5 ergs cm -2 
sec -1 Hz -1 sr _1 over the redshift range 1.7 < z < 3.8. Sev- 
eral other authors carried out this analysis on other data 



sets (iLu et alj 



Willi ger et al.l 



1996) 



Lu et al. 



1991 



1994 



Kulkarni fc Fall 



19931 ; IBechtoldl 1 1994 



Cristiani et al.l 19951 ; iGiallongo et al 



(1991) found the same result of lBaitlik et al 



( 1988 ) in the s ame r edshift interval but using 38 QSO spec- 
tra. IBechtoldl (|l994), with 34 low resolution QSO spectra 
covering the range 1.6 < z < 4.1, fou nd log Jll — — 20.5 . 
Thi s value, 3 times larger than that of lBaitlik et alj (Il988l) 
andlLu et all l|l99ll ). was not confirmed by IGiallongo et al.l 
jl996h They obtained log Jll ^ -21.3^;^, using 10 higher 
resolutio n (R ~ 25000) Q SO spectra in t he same redshift in - 
terval as IBechtoldl l|l994l ). More recently. IScott et al.l (l200d ) 
considered a sample of 74 intermediate resolution QSO spec- 
tra in the redshift interval 1.7 < z < 4.1, from which they 
obtained a value log Jll — — 21.11q'27. These authors paid a 
particular attention to the correct estimate of the systemic 
redshift of the QSOs. Indeed, if a QSO redshift lower than 
the true one is used in the analysis, the ionising effect of the 
QSO on the Lya clouds, and thus the derived value of Jll, 
are over-estimated. 

In the standard analysis of the proximity effect it is 
assumed that the matter distribution is not altered by the 
presence of the QSO. The only difference between the gas 
close to and far away from the QSO is the increased pho- 
toionisation rate due to the QSO emission. A consequence 
of this hypothesis is that there should be a correlation be- 
tween the strength of the proximity effect and the luminos- 
ity of the QSO. However, ob servational results are not con- 
clusive on this subject (e.g. ILu et al.lll99ll; iBechtold 



. , _ - - -jll — 

ISrianand fc Kharell 996: but see also lLiske fc" Willigerll200l 
It is in fact likely that QSOs occupy over-dense regions. Hi- 
erarchical models of structure formation predict that super- 
massive black ho les, that are though t to power QSOs, are 
in massive halos dGranato et al. | |2004 iFontanot et al. 2006; 
Ida Angela et al.l |200S| ) which are strongly biased to high- 
density regions. 

The main aim of this work is to investigate the density 
distribution of matter close to QSOs from the proximity 
effect. To this purpose, we applied the FLO (From L i nes to 
Over-densities) technique, developed in ISaitta et al. I l|2007l . 
Paper I), to a sample of 21 high resolution, high signal-to- 
noise ratio QSO spectra. FLO converts the list of H I column 
densities of Lya lines in a QSO spectrum into the underlying 
tih hydrogen density field. This method significantly reduces 
the drawbacks of the line fitting approach, in particular, the 
dependence on the fitting tool, the subjectivity of the result, 
and the strong dependence of statistics on the number of 
weak lines, which is in general poorly known. We could also 
constrain the value of the UVB ionisation rate, TubvEI, by 
matching the recovered density field in the Lya forest region 
with the mean cosmic density. 

Similar analyses were performed in previous works us- 



1994; 



1 The Lyman limit corresponds to the hydrogen ionisation energy 
E ion = 13.6 eV or A LL = 912 A. 

2 Fuvb = 4-7T / a v {J v lhv)dv where, <t„ is the photo- 
ionisation cross section, i/ll is the frequency at the Lyman limit 
and J„ = JllWz'll) -7 - 



ing a different approach, based on the determination of the 
cumul ative probability distri bution function of pixel optical 
depth. iRollinde et alj l|2005h studied the density structure 
around QSOs using 12 high resolution-spectra that belong 
also to our sample. These authors marginally detected the 
presence of an over-density at separations 3 r ^ 15 proper 
Mpc, assuming an hydrogen ionisation ra te Tuvb = 10~ 12 
(corresponding to log Jll — —21.4). iGuimaraes et al] 



i|2007l ). using the same technique, investigated the distribu- 
tion of matter density close to 45 high-redshift (z cm ~ 3.8) 
QSOs observed at medium spectral resolution. Their study 
reveals gaseous over-densities on scales as large as ~ 15 Mpc, 
with higher over-densities for brighter QSOs. 

The paper is organised as follows: in Section [2j we in- 
troduce the FLO technique and describe an improvement to 
the method. Section [3] reports the details of the observed 
data sample, the QSO emission redshift and luminosity de- 
termination, and the fitting analysis. In SectiorJJ] the prox- 
imity effect of QSOs on the surrounding gaseous medium 
is pointed out and the over-density close to QSOs is recon- 
structed with FLO. The characteristics of the cosmological 
simulations used to obtain the mock spectra and the com- 
parison with observations are described in Section 6. Our 
conclusions and the prospect for the future are listed in Sec- 
tion 7. 

Throughout this paper, we adopt a ACDM cosmol- 
ogy with the following values for the cosmological param- 
eters at 2 = 0: Q m = 0.26, fi A = 0.74, Q b = 0.0463, 
n 3 = 0.95, a 8 = 0.85 and Ho = 72 km s" 1 Mpc" 1 . These 
parameters are consistent with the best fits values obtained 
from the latest results on the cosmic microwave background 
l|Komatsu et alj 20081) and from other flux statistics of the 
Lya forest (e.g. lViel et"aill2006l) . 



2 THE FLO TECHNIQUE 

FLO (From Lines to Over-densities) was introduced in Paper 
I. The physical hypotheses at the base of this procedure are 
briefly described in the following paragraphs. 

Traditionally, the analysis of the Lya forest was based 
on the identification and Voigt fit of the absorption lines in 
order to derive the central redshift, the column density and 
the Doppler parameter (measuring the velocity dispersion 
in the line). This approach has two main drawbacks: (i) the 
subjectivity of the decomposition into components: the same 
complex absorption can be resolved by different scientists 
(or software tools) in different ways, both in the number 
of components, and in the values of the output parameters 
for a single component; (ii) the blanketing effect of weak 
lines: they can be hidden by the stronger lines, so that their 
exact number density is unknown and has to be inferred 
from statistical arguments. Unfortunately, since the weak 
lines are also the most numerous, the uncertainty in their 
exact number is transformed into a systematic error of the 
computed statistical quantities. 

The FLO technique extends the line fitting approach by 
identifying a new statistical estimator describing the phys- 
ical properties of the underlying IGM, the hydrogen den- 
sity riH , which is link ed to the mea sured H I column density 
through the formula (Schaye 2001): 



AT (HI) 



3.7 x 10 13 



(1+5) 1 
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where, 5 = nu/inn) — 1 is the density contrast, To, 4 = 
To/10 4 K is the temperature at the mean density, Ti2 = 
Tuvb/IO -12 s _1 is the H photo-ionisation rate due to the 
UV background, f g ~ f2b/^m is the fraction of the mass 
in gas and a is the index of the temperature-density rela- 
tion for the IGM which depends on the ionisation history of 
the Universe. Equation [1] relies on three main hypotheses: 
(i) Lya absorbers are close to local hydrostatic equilibrium, 
i.e. their characteristic size will be typically of the order of 
the local Jeans length (Lj), which can be appr oximated a^f] 
jNusser fc Haehneltllioool ; IZaroubi et al.ll2006l ): 

Lj (2) 

(^y 1/2 (%t) 1/2 my /2 mr /2 

in comoving Mpc, where h = f/o/100 km s _1 Mpc -1 and 
the other parameters were already defined; (ii) the gas is in 
photo- ionisation equilibrium ; (iii) the 'effective equation of 
state' (|Hui fc Gnedinlll997r i. 

T = T (5 + l) a , (3) 

holds for the optically-thin IGM gas. 

In order to apply eq.[T]we have, first of all, to go through 
the Voigt fitting process of the Lya forest absorptions in a 
QSO spectrum. Then, to transform the list of Hi column 
densities of Lya lines into the matter density field which 
generated them, we have to perform the following steps: 

(1) group adjacent Lya lines into absorbers of size of 1 Lj 
with column density equal to the sum of column densities 
and redshift equal to the weighted average of redshifts, using 
column densities as weights. The absorbers are created with 
a friend-of-friend algorithm: 

(i) the spatial separation between all the possible line 
pairs is computed and the minimum separation is compared 
with Lj, computed at the N(R l)-weighted redshift mean of 
the pair; 

(ii) if the two lines of the pair are more distant than the 
local Lj, they are classified as two different absorbers, stored 
and deleted from the line list; 

(iii) if the two lines are closer than the local Lj, they are 
replaced in the line list by one line with a redshift equal to 
the iV(H l)-weighted mean of the two redshifts and a column 
density equal to the sum of the two column densities; 

(iv) the procedure is iterated until all the lines are con- 
verted into absorbers. 

(2) transform the list of column densities of absorbers into 
a list of 5 inverting eq. [T] 

(3) bin the redshift range covered by the Lya forest into 
steps of 1 Lj and distribute the absorbers onto this grid, 
proportionally to the superposition between absorber size 
(which is again 1 Lj) and bin. Two cases are considered 
for the treatment of the empty bins: in the 'lower limit' 
case the bin is filled with an absorber of null column density 
(corresponding to 5 — —1), while in the 'upper limit' case the 
bin is filled with an absorber with hydrogen density contrast 

3 Note that the formula in the original papers has the wrong sign 
for the exponent of a + 1. 




log(<5+ I) 

Figure 1. Distribution of the fraction of bins as a function of their 
log(5 + l)-value for the true field (solid line) and the recovered 
one (dotted line) from the simulated lines of sight in the upper 
limit case. The smoothing scale is the Jeans length. 
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Figure 2. Contour scatter plot of the true versus reconstructed 
<5 fields from simulations in the upper limit case. The smoothing 
scale is the Jeans length. The contours show the number density 
of pixels which increases by a factor of 10 at each level. 



corresponding to the minimum detectable column density in 
our data, logiV(Hl) = 12 cm -2 , at the redshift of the bin. 

For the IGM parameters in eqs.[T]and[2]we adopted the 
same values of paper I: To, 4 = 1.8, Ti2 = 1 and a = 0.6. 

2.1 An improvement of the FLO method 

In paper I, we tested the quality of the density reconstruc- 
tion of FLO using synthetic QSO spectra drawn from a 
cosmological simulation (the same that is described in Sec- 
tion 6). Figure [T] and [2] shows the results of those tests. Both 
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Figure 3. As in Fig. [T] but using a filtering scale Lp ~ 0.7 Lj. 



the true density field and the reconstructed one along the 
simulated lines of sight were rebinned into steps of 1 Jeans 
length . The recovered density field from the simulated lines 
of sight presents two main problems: (i) most of the points 
of the true density field belonging to under-dense regions are 
not recovered with the correct <5-value, they are accounted 
for as empty bins or, for those approaching the average den- 
sity, their value is over-estimated and they are moved to 
moderately over-dense regions; (ii) the number of points in 
the over-dense regions is over-estimated due to a systematic 
assignment of larger-than-true (5-values to those points. The 
latter effect causes also an over-estimate of the average 8 of 
the distribution which is 40 — 50 percent higher than that 
of the true density field. In paper I, we solved this problem 
by normalising both the true and the reconstructed <5-field 
in order to have the same mean value, {5 + 1} = 1.0. This 
solution is not viable for the present analysis, so we looked 
for the physical reasons of the wrong reconstruction in order 
to improve the performances of FLO. 

The failure in reproducing the under-dense regions is 
likely due to the fact that gas below the average density is 
still expanding and our primary hypothesis, the local hy- 
drostatic equilibrium, cannot be applied. There is no simple 
solution to this problem, however we have already tested in 
paper I that under-dense regions have a negligible effect on 
statistical quantities. 

The problem regarding the moderately over-dense re- 
gions depends on the scale adopted for the reconstruction 
of the absorbers, the Jeans length Lj, which appears to be 
too large. As a consequence, absorbers have on average too 
large column densities which translate into an over-estimate 
of the 5-values for the over-dense regions. The need for a 
smoothing length smaller than Lj to better reproduce the 
density field traced by the L yq forest, was discussed in de- 
tail by iGnedin fc Huil (|l998h . They affirmed that the cor- 
rect filtering scale depends on the re-ionisation history of 
the Universe and is in general smaller than the Jeans scale 




iog(5 true +i) 

Figure 4. As in Fig. [2] but using a filtering scale Lp ~ 0.7 Lj. 
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Figure 5. Lyc? forest redshift coverage of the QSOs in our sample. 



after re-ionisat i on, a nd larger prior to it. We used eq. A4 of 
IGnedin fc Huil (Il998l ) to compute Lp, adopting z r ci = 11 and 
z bs = 2.5. The results using the filtering scale Lf ~ 0.7 Lj 
are shown in Fig. [3] and [4] The moderately over-dense re- 
gions are now successfully recovered both in the number of 
points and in the 5-values, and the average 5 of the distribu- 
tion is in agreement with the true one within 10 percent. We 
verified that varying the re-ionisation redshift or the temper- 
ature at the mean density by ~ 20 percent does not have a 
significant effect on the density reconstruction process. 

In the following analysis, the filtering scale Lf replaced 
the Jeans length in the FLO algorithm described before. 
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Table 1. Relevant properties of the QSOs forming our sample. See text for further details. 



QSO 


z em (Ref.) 


line 


Az Lya 


bj 


r F 


F07 


r/io 40 

T02 


HM 


log ^LL 


(Mpc) 


HE1341-1020 a 


2.142(1) 


Mgn 


1.6599-2.142 


18.68 


17.52 


0.11±0.03 


0.10±0.02 


0.12±0.02 


30.712 


2.9 


Q0122-380 


2.2004(2) 


113 


1.709-2.2004 


17.34 


16.70 


0.8±0.3 


0.5±0.1 


0.6±0.1 


31.267 


5.5 


PKS1448-232 


2.224(1) 


Mgn 


1.729-2.224 


17.09 


16.87 


1.0±0.3 


0.56±0.09 


0.6±0.1 


31.377 


6.3 


PKS0237-23 


2.233(1) 


Mg n 


1.737-2.233 


16.61 


16.21 


1.3±0.6 


0.8±0.2 


1.0±0.2 


31.573 


7.9 


J2233-606 


2.248(1) 


1 


1.7496-2.248 


16.97 


17.01 


1.4±0.2 


0.66±0.09 


0.8±0.1 


31.437 


6.7 


HE0001-2340 


2.265(1) 


Mgn 


1.764-2.265 


16.74 


16.46 


1.2±0.5 


0.7±0.2 


0.9±0.2 


31.538 


7.6 


HE1122-1648'' 


2.40(3) 




1.878-2.344 


16.61 


16.32 


1.8±0.5 


1.0±0.2 


1.2±0.2 


31.679 


8.9 


Q0109-3518 


2.4057(4) 


[O m] 


1.883-2.4057 


16.72 


16.37 


1.8±0.5 


1.0±0.2 


1.1±0.2 


31.640 


8.5 


rlrjZzl t-Zoio 


Z.41z^l I 




l.ooo-z.ooo 


ID. 4 / 


i ft i p. 


Z. 1 ±U.D 


l.o±U.Z 


1 KJ-(\ o 

l.D±U.z 


Q 1 7/1/1 

o 1 . / 44 


O ft 


Q0329-385 


2.437(5) 


Mgn 


1.9096-2.437 


17.20 


16.91 


1.6±0.3 


0.7±0.1 


0.8±0.1 


31.472 


7.0 


HE1158-1843 a ' 6 


2.448(1) 




1.919-2.391 


17.09 


16.84 


1.3±0.2 


0.70±0.09 


0.8±0.1 


31.525 


7.4 


HE1347-2457 a 


2.5986(4) 


H/3 


2.046-2.5986 


17.35 


16.14 


0.7±0.3 


0.7±0.1 


0.8±0.1 


31.542 


7.6 


Q0453-423 


2.669(1) 


Oi 


2.106-2.669 


17.69 


16.74 


0.7±0.4 


0.5±0.1 


0.6±0.1 


31.451 


6.8 


PKS0329-255 6 


2.696(1) 




2.129-2.635 


17.88 


17.71 


0.8±0.1 


0.38±0.06 


0.44±0.07 


31.396 


6.4 


HE0151-4326 6 


2.763(1) 




2.186-2.701 


17.48 


16.93 


1.7±0.4 


0.9±0.1 


1.0±0.1 


31.594 


8.1 


Q0002-422 


2.769(1) 


Oi 


2.191-2.769 


17.50 


16.89 


1.5±0.4 


0.8±0.1 


1.0±0.1 


31.589 


8.0 


HE2347-4342 a 


2.880(1) 


Oi 


2.285-2.880 


17.12 


16.30 


2.0±0.7 


1.3±0.2 


1.5±0.2 


31.810 


10.4 


HS1946+7658 a 


3.058(6) 


Oi 


2.435-3.058 


16.64 


15.76 


5±1 


2.7±0.4 


3.1±0.5 


32.118 


14.8 


HE0940-1050 


3.0932(4) 


H/3 


2.465-3.0932 


17.08 


16.08 


3±1 


2.0±0.3 


2.3±0.3 


31.972 


12.5 


Q0420-388 a 


3.1257(1) 


Oi 


2.493-3.1257 


17.44 


16.70 


2.6±0.5 


1.2±0.2 


1.4±0.2 


31.848 


10.8 


PKS2126-158 


3.292(7) 


[Om] 


2.633-3.292 


17.54 


16.37 


4±1 


1.8±0.3 


2.0±0.3 


31.958 


12.3 



a QSOs with associated absor ption systems; b QSO no t considered in t he proximity study. 

References: (1) this pa per: (21 [Sulentic et al.ll2004l; (31 iKim et"alll2002l : (4) P. Marziani, private comm.; (5l lEspev et al.lll989l : (6) 
iFan fc Tvtlenll994l : (71 IScott et alll2000l . 



3 OBSERVED DATA SAMPLE 

Most of the observational data used in this work were ob- 
tained with the UVES spectrograph |Dekker et al.ll200ol ) at 
the Kueyen unit of the ESO VLT (Cerro Paranal, Chile) 
in the framework of the ESO L arge Programme (LP ): "The 
Cosmic Evolution of the IGM" (jBergeron et al.ll2004h . Spec- 
tra of 18 QSOs were obtained in service mode with the aim 
of studying the physics of the IGM in the redshift range 
1.7-3.5. The spectra have a resolution R ~ 45000 and a typ- 
ical signal to noise ratio (SNR) of ~ 35 and 70 per pixel at 
3500 and 6000 A, respectively. The wavelength range goes 
from 3000 to 10,000 A, except for two intervals of about 100 
A, centred at ~ 5800 and 8600 A where the signal is ab- 
sent, due to the gap between the two CCDs forming the red 
mosaic. In the spectra of the two QSOs at higher redshift 
(Q0420-388 and PKS216-158) there are three gaps centred 
at ~ 5640 and 8600 A of about 100 A and at ~ 6680 A of 
wi dth ~ 50 A. Detail s of t he data reduction c an be found 
in IChand et al.l (|2004l ) and lAracil et alj l|2004l ). In particu- 
lar, the continuum level was estimated with an automatic 
iterative procedure which underestimates the true contin- 
uum in the Lya forest to about 2 percent at z ~ 2.3. In 
the process of fitting the lines in the Lya forest, we cor- 
rected the continuum level in the spectral intervals where it 
was clearly underestimated by interpolating the regions free 
from absorption with polynomials of 3rd order. 

We added to the main sample 3 more QSO spectra with 
comparable resolution and SNR: 

- J2233-606 (|Cristiani fc D'Odoricol l200Ch . Data for this 
QSO were acquired during the commissioning of UVES in 
October 1999. 

- HE1122-1648 (|Kim et al.ll2002l ). Data for this QSO were 



acquired during the science verification of UVES in February 
2000. The reduced and fitted spectrum was kindly provided 

to us by Tae-Sun Kim. 

- HS1946+7658 (|Kirkman fc Tvtlerl li"997T ). Data for this 
QSO were acquired with Keck/HIRES in July 1994. 

In the following sub-sections, the procedure to derive 
the QSO emission redshift and luminosity and the Lya line 
lists is described. Table Q] reports for each QSO in the sam- 
ple the final emission redshift and how it was computed, 
the studied Lya redshift range, the apparent magnitude of 
the QSO, and the corresponding H ionisation rate resulting 
from the three adopted QSO spectra, the luminosity at the 
Lyman limit and the radius of influence of the QSO ionis- 
ing flux. Figure shows the distribution in redshift of the 
Lya forests for all the QSOs of the sample. We considered 
the range between 1000 km s^ 1 red-ward the Ly/3 emission, 
to avoid contamination by Ly/3 absorption lines, and the 
Lya emission. 



3.1 Estimate of the QSO systemic redshifts 

The knowledge of the correct systemic redshift of the QSO is 
of fundamental importance when using the proximity effect 
to estimate both the intensity of the UV ionising background 
and the density structure close to the QSO itself. 

Emission redshift of QSOs at z Gm £ 1.5 are gener- 
ally computed by the positions of the most prominent 
UV emission lines, in particular Hi Lya A 1216 and 
C iv A 15 49. However, i t was assessed by several stud- 
ies (e.j 

"reenl 



Gaskelll [1983: IWilkesI [l986l: lEspev et al.l Il989l ; 

19921; 



Corbinfll990l; iTytler fc Fan! Il992l; iBoroson fc Gree n tM>> 
Laor et al.ll 19951 ; iMarziani et al.lll996l ; lMcIntosh et al.|[l99 
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IVanden Berk et alj l200ll ; ISulentic et ail 120041 ) that high- 
ionisation emission lines (e.g., C IV, N v A 1240, C ill] A 1909, 
and Hi Lya ) are on average blue-shifted by several hun- 
dreds of km s -1 with respect to low-ionisation lines (e.g., 
Ol A1304, Mgn A2798) and the permitted Hi Balmer se- 
ries. In contrast, redshifts from narrow forbidden lines (e.g., 
[O 11] A 3727, [O in] A 5007) are observed to be within 100 
km s -1 of the broad Mgn and Balmer lines. Furthermore, 
in local active galactic nuclei, redshifts from narrow forbid- 
den lines showed agreement to ~ 100 km s _1 of the ac- 
cepted systemic frame determined by stellar absorp tion fea- 
tures and Hi 21 cm emission i n the host galaxies |Gaskelll 
1 19821 ; IVrtilek fc Carletonlll985l ; iHutchings et al.lll987l) . 

In order to estimate the correct redshift, we carried out 
a detailed analysis of the emission lines of the QSOs in our 
sample both using data from the literature and directly fit- 
ting the lines in the UVES spectra. In decreasing order of 
precision, we adopted the redshifts measured by: i) narrow 
forbidden lines, mainly [Oiii], ii) H/3, iii) Mgn, and iv) Ol. 
There are 5 QSOs in the sample for which none of these lines 
was measured. We excluded those objects from the proxim- 
ity study, disregarding the portion of the spectrum within 
5000 km s -1 of the best estimate of the emission redshift 
(generally obtained from Lya and/or Civ emission lines). 

The determination of the redshifts from the UVES spec- 
tra were obtained by re-binning the region of the emission, 
normalising it to the local continuum and fitting a Gaussian 
profile to the line. Table [T] gives the QSO redshifts and the 
details of the estimate. 

3.2 Estimate of the QSO ionising fluxes 

Getting closer and closer to the QSO, the UV ionising field 
becomes dominated by the intrinsic QSO emission flux. In 
order to derive the matter density distribution around the 
QSO using the observed variation of the absorption features 
in the QSO spectrum, a reliable determination of the intrin- 
sic luminosity of the QSO has to be obtained. 

Magnitudes of the obj ects in our sample were taken 
from the GSC-II catalogue (|McLean et al.ll2000h and are re- 
ported in Tab.[T] In order to estimate the corresponding in- 
trinsic bolometric luminosity, we adopted the following pro- 
cedure. 

We consid e red t he QSO template library defined in 
iFontanot et all (|2007l ). based on high quality SDSS QSOs 
spectra in the redshift interval 2.2 < z < 2.25. This redshift 
range was chosen in order to maximise the level of complete- 
ness of the sample and the wavelength interval long-wards 
of the Lya emission. Moreover, in this redshift range the dy- 
namical response of the SDSS spectrograph is such that the 
Lya line is completely sampled in all spectra. In the original 
paper, the authors considered the rest-frame spectra of the 
215 QSOs forming the final sample and they used a con- 
tinuum fitting technique in order to extend the information 
blue- ward of the Lya . A mean continuum slope 7 = 0.7±0.3 
was obtained for the objects in the library, where /„ oc v~"< 
and /„ (corresponding to An J v ) is the QSO flux in units of 
ergs cm -2 s -1 Hz -1 . IFontanot et al.l (|2007l ) demonstrated 
that this library is suitable for predicting QSO colours up 
toz~ 5.2. 

We used the template spectra in the library to com- 
pute a synthetic bj and tf magnitude at each emission 



redshift listed in Table [T] For reproducing the bj and tf 
photograp hic magnitudes, the response of t he spectral code 
PEGASE (|Fioc fc Rocca-Volmerangel [l997) was assumed. 
Then, the templates were renormalised in each band sep- 
arately, by requiring the synthetic magnitude to match the 
observed one, and the renormalised spectra were used to 
give a prediction for the AB magnitudes at 912 A {M\{ 2 an d 
Mgf 2 respectively) . The quantity AM912 = M§( 2 ~ ^912 was 
adopted as an estimator of the agreement between the slope 
of the template and the intrinsic slope of the considered 
QSO: we then associated to each observed QSO the template 
with the smaller AM912. For 18 out of 21 QSOs in our sam- 
ple this procedure gave AM912 values lower than 0.001 mag. 
For the remaining 3 objects the values were respectively 0.1 
(PKS0329-255), 0.2 (HE1347-2457) and 0.4 mag (HE1341- 
1020). In these cases, there was no template in the library 
which reproduced the correct intrinsic slope of the observed 
QSO. The AM912 can be taken as a measure of the system- 
atic error on the luminosity computed for these 3 objects. 

The selected and renormalised template spectra were 
completed in the region blue- ward of the Lya emission using 
three possible extrapolations: 

F07: the continuum slope of the template spectrum red- 
ward of the Lya emission; 

T02 : a fixed power la w with slope 7 = 1.8 ( |Madau et al.l 
1 19991 : iTelfer et al.ll2002j ) ; 

HM: a fixed power law with slope 7 = 1.5 
l|Haardt fc Madaulll99rJ) . 

Then, the spectra were used to estimate the QSO monochro- 
matic luminosities, L v (ergs s _1 Hz -1 ). The hydrogen 
photo-ionisation rates due to the QSO radiation were ob- 
tained through the formula: 

TqsoIs- 1 ] = (4) 

where r a b s is the distance (in cm) between the QSO emis- 
sion redshift and the absorber and F is obtained by integrat- 
ing the monochromatic luminosity in the wavelength range 
300 < A < 912 A: 

r= / ^a v dv (5) 

/ hv 

where a v is the ab sorbing cross-section for neutral hydrogen 
i|Osterbrocklll989l ). 

The error on F was computed by taking into account 
the uncertainties in the QSO magnitudes. We randomly ex- 
tracted 100 values of bj and tf in the corresponding al- 
lowed range and we repeated the previous procedure. The 
obtained mean F values are in good agreement with the es- 
timate based on the observed magnitudes, and we adopted 
the variance on the 100 realizations as the error on V. The 
values of F with their uncertainties, for the three investi- 
gated QSO continua, are reported in Table [1] The three 
estimates F07, T02 and HM agree within 2a. In particular, 
the r for the T02 and HM models are always within 1 a. As 
a consequence the differences between the results obtained 
adopting the three Y values could be interpreted also as due 
to the uncertainties in the F itself. 

The radius of the sphere of influence of each QSO in 
our sample is obtained by relating the intensity of the UV 
ionising background at the Lyman limit, Jll, and the lumi- 
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nosity of the QSO at the same frequency, Lll, measured as 
described above, 



4^(Lll/Jll 



,V2 



(6) 



The luminosities depend slightly on the adopted slope, how- 
ever the differences are small so we used the average of the 
three values to compute r eq . For the UV background, we 
adopted the value Jll ~ 4 x 10 -22 ergs cm -2 s _1 Hz -1 
sr~\ corresponding to an ionisation rate Fuvb — 10 -12 s _1 . 
The resulting average Lll and the radii expressed in proper 
Mpc are reported in Table [1] 



3.3 Compilation of the line lists 

All the lines in the Lya regions of the LP QSOs 
plus J2233-606 were fitt ed with the FITLYMAN tool 
|Fontana fc Ballestejfl995h of the ESO MIDAS data reduc- 
tion packag43- In the case of complex saturated lines, we 
used the minimum number of components to reach \ 2 ^ 1-5. 
Whenever possible, the other lines in the Lyman series were 
used to constrain the fit. The minimum Hi column density 
detectable at 3 a with the SNR of the spectra of our sample 
is logiV(Hl) ~ 12 cm -2 . 

Metals in the forest were identified and the correspond- 
ing spectral regions were masked to avoid effects of line blan- 
keting. The same treatment was given to Lya lines with 
column density AT (Hi) ^ 10 16 cm -2 since, on the one hand, 
strong H I lines can hide weaker lines as much as metal lines, 
and on the other hand the application of the FLO algorithm 
is valid only in the linear or mildly non-linear regime, or for 
over-densities 5 ~ few x 10. We eliminated Lya lines with 
Doppler parameters b ^ 10 km s~ , that are likely unidenti- 
fied metal absorptions. They represent about the 7 percent 
of the total sample of Lya lines. The output of this analysis 
is a list of Lya lines for each QSO with central redshift, H I 
column density, Doppler parameter and the corresponding 
errors obtained with FITLYMAN. 

The Lya forests of the remaining two QSOs, HE1122- 
1648 and HS1946+7658, were fitted with the VPFIT0 pack- 
age. The same constraints on Doppler parameter and column 
density were applied to obtain the final Lya line lists for 
these QSOs. The difference in the FLO results due to the dif- 
ferent fitting software, FITLYMAN and VPFIT, is negligible 
and was discussed in Paper I. 



4 THE PROXIMITY EFFECT 

Moving closer to the QSO a decrease of the number density 
of Lya absorption lines with respect to the mean value is 
expected. 

We applied FLO to recover the density field in the neigh- 
bourhood of the QSOs in our sampl^fl, deliberately neglect- 
ing the contribution of the QSO radiation field to the UV 
ionising flux. The density field traced by the Lya forest was 
computed in each QSO spectrum of our sample. The data 



4 http://www.eso.org/midas 

5 http:/ /www. ast. cam. ac.uk/~rfc/vpfit. html 

6 Note that due to the uncertainty in the emission redshifts, 5 
QSOs were excluded from the computation of the proximity effect. 
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Figure 6. The mean of the S field reconstructed around the QSOs 
in our sample, neglecting the effect of the UV flux from the QSO 
themselves. The solid dots represent the upper limit case while 
the open stars trace the lower limit case (see text). The shaded 
area marks the 3 a region for the upper limit case, the dotted lines 
mark the corresponding 1 and 2 a levels and the dashed line the 
average IGM density contrast. The lower panel shows the number 
of objects contributing to each bin. 



were re-binned into bins of 8 proper Mpc in length. Bins 
covered by more than 33 percent by masked intervals were 
eliminated from the final count for the single object. Then, 
in each bin, the mean of all the S values contributing to 
that bin (one per QSO at maximum) was computed. The 
resulting 5-field is shown in Fig. [S] together with the num- 
ber of QSOs contributing to each bin. The upper and lower 
limit case are shown, the difference between the two recon- 
structions is negligible. The dashed line in the upper panel 
represents the average of S values in the upper-limit case at 
separations larger than 20 Mpc, no longer affected by the 
ionising flux from the QSO. The dotted lines are the 1 and 
2 a standard deviations and the shaded region represents 
the 3 a level. A significant decrease of the density field is 
observed in the first bin at proper separation from the emit- 
ting QSO of less than 8 Mpc, corresponding to the average 
radius of influence of the QSOs in our sample. 



4.1 The density structure around QSOs 

In order to recover the density field in the region close to 
the emitting QSO, its ionising flux has to be added to the 
UV background, and the total ionisation rate has to be used 
in eq. [T] of the FLO algorithm. We define 



r' 12 = 



Tuvb + FqsQ 

io- 12 



(7) 



where Fqso was defined in eqj4] As explained in Section ^. 21 
the value of Tqso depends on the adopted slope for the 
QSO continuum in the region blue- ward of the Lyman limit. 
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Figure 7. Mean 8 field in the upper (top panels) and lower (bottom panels) limit case in the region close to the emitting QSO (which 
is located at r = 0.0) for the QSO spectra in our sample. In each row, the three panels differ for the adopted QSO continuum slope (see 
Section l3.2l l. The three curves in each panel differ for the adopted value of Ti2. The 3<r fluctuations in the IGM region (far away from 
the QSO) are represented by the shaded area. 



We investigated how the three different assumptions for the 
slope influenced the final result. 

In Fig. [7J we present the mean 8 field for the lower 
and the upper limit case recovered from the spectra in our 
sample in the region within 20 Mpc from the radiating source 
in bins of 4 proper Mpc. We studied how the FLO output 
was modified by assuming different ionisation rates for the 
UV background, Ti2 = 1, 2, 3, and different slopes for the 
QSO continuum. 

The first important result of our analysis is that FLO 
recovers the over-density hosting the QSOs in each one of 
the panels of Fig. [7] The second result is that the value of 5 
for the peak varies significantly depending on several factors 
that are discussed in the following. 

1. Treatment of the empty bins: at variance with the region 
away from the QSO, the lower and upper limit case in the 
treatment of the empty bins makes a difference of a factor 
of ~ 2 for the 5 value of the peak of the over-density where 
the QSO resides. This is due to the QSO strong ionising 
field that, when applied in eq. [TJ transforms even small col- 
umn densities into large over-densities. Increasing the SNR 
of our spectra would decrease the minimum detectable H I 
column density, thus decreasing the difference between the 
two limiting cases. 

2. QSO continuum slopes: the difference between the two 



extrapolations with the standard fixed power laws, with in- 
dexes 7 = 1.5 and 1.8, is ~ 10 percent, while there is an 
increase in the peak value between ~ 30 and 45 percent 
when the F07 slopes are adopted. In the latter model, the 
continuum in the blue is obtained by extrapolating the slope 
red-ward of the Lya emission, so it can be considered as a 
sort of upper-limit to the true continuum, that cannot be 
recovered due to contamination by the IGM absorption. 

3. Ionisation rate of the UVB: this parameter does not 
have any influence on the S value for the QSO peak, as ex- 
pected. On the other hand, the effect on the average value 
of the (5-field away from the QSO is significant and will be 
discussed in Section H. 3 1 

4. QSO systemic redshifts: the uncertainties in the emis- 
sion redshifts of the QSOs whose spectrum was used to de- 
termine the over-density are of the order of 100 — 200 kms -1 
corresponding to ~ 400 — 800 proper kpc at z ~ 2.5. As 
a consequence, we do not expect a major influence on the 
value of the QSO over-density which is computed in a bin 
of 4 proper Mpc. A significant improvement in the result 
would be obtained if the systemic redshift of all the QSOs 
in the sample would be measured with intermediate resolu- 
tion spectra in the infrared. 

5. Associated absorption systems: metal absorption sys- 
tems within ~ 5000 km s" 1 (corresponding to ~ 20 proper 
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Figure 8. Comparison between the recovered average density 
field for the total sample (solid dots) and the sample without 
the QSOs with associated systems (empty stars) in the reference 
model with the F07 slopes. The shaded regions mark the 3 <r inter- 
val an the dotted lines the 1 and 2 a level for the reduced sample. 
In the lower panel the number of objects for the two samples is 
reported. 



Mpc) from the emission redshift of the QSO are defined as 
associated to the QSO. Some of these systems could indeed 
be intrinsic to the QSO, that is they could be very close to 
the emitting region and their position could be determined 
only by the large velocity at which they were ejected from 
the QSO itself. Counting these lines as intervening would 
result in an over-estimate of the density field in the prox- 
imity of the QSO. Six QSOs in our sample show associated 
absorption lines (AAL). Figure [5]reports the comparison be- 
tween the recovered density field of the total QSO sample 
and of the sample without AAL for the reference model and 
the F07 slopes. There is a decrease of ~ 40 and 15 percent in 
the 8 value of the peak in the lower and upper limit-case, re- 
spectively. A slight increase in the scatter of the IGM density 
field is also measured due to the lower statistics. The point 
at a separation of 98 Mpc is above the 3 a level, however 
the fluctuation is due to a single QSO, PKS1448-232 and 
disappears when the QSO is eliminated from the sample. 



4.2 Comparison with previous results 

The density distribution in the region close to bright QSOs 
was studied in two p revious works (|Rollinde et all [2005; 
ICuimaraes et"alll2007h with a different technique: the vari- 
ation of the cumulative probability distribution function 
(CPDF) of the optical depth measured in the Lya forest 
region away from and close to the QSO. The advantage of 
this method with respect to FLO is that the fitting of the 
lines is not necessary since the observed quantity used is 
the pixel-by-pixel transmitted flux. On the other hand, be- 
fore investigating the change of the CPDF with the distance 
from the QSO, it is necessary to evaluate its evolution with 



redshift and subtract it. Then, the change in the CPDF has 
to be translated into a va riation of the densit y field with 
a bootstrap technique (see iRollinde et al. 2005,' for the de- 
t ail ed description). 

IRollinde et all l|2005l ). using a subsample of our sam- 
ple formed by 12 QSO spectra with (z e m) — 2.5, found a 
significant over-density for separations between ~ 3 and 15 
proper Mpc assuming Fi2 < 3. They assumed in the com- 
putation a QSO continuum with a fixed power law slope 
7 = 0.5 and a temperat ure-density relation index a = 0.5. 
ICuimaraes et all l|2007h applied the same technique to a 
sample of 45 QSO spectra at intermediate resolution and 
at redshift (z em ) ~ 3.8. An over-density extending to sep- 
arations of ~ 15 proper Mpc is detected adopting the pa- 
rameters, Ti2 = 1, To. 4 = 1, a = 0.0 and 7 = 0.7. These 
authors claimed also the detection of a correlation between 
the over-density and the luminosity: brighter QSOs reside 
in higher over-densities. 

The main difference between our result and the previ- 
ous ones is that we clearly detected an over-density limited 
to a few Mpc from the QSO, while the results of the optical 
depth analysis showed a smooth decrease from the peak to 
the IGM density extending for more than 10 Mpc. In par- 
ticular, the difference w ith the result at higher redshift by 
ICuimaraes et all (|2007T l cannot be ascribed to a difference 
in QSO intrinsic luminosity, since the average luminosity of 
the two samples is the same: logl/LL — 31.7. We will see 
in Section [5] that the cosmological hydro- simulations that 
we used for comparison show very narrow density peaks as 
found in our recovered density field. Very s imilar results were 
obtained bv lFaucher-Giguere et all (|2008l ) who investigated 
with hydrodynamical simulations the bias introduced by the 
QSO over-density in the estimate of the UV background 
intensity from the proximity effect. They considered three 
ranges of masses 3.0 ± 1.6 x 10 11 M Q h' 1 , 3.0 ± 1.6 x 10 11 
M h' 1 and 6.0 ± 3.2 x 10 12 M h' 1 for the DM haloes 
hosting QSOs in the simulations and averaging the results 
of 100 lines of sight obtained overdensity profiles extending 
to ~ 3 — 5 proper Mpc at redshifts z ~ 2 — 3. 

4.3 Constraints on the IGM physical parameters 

The FLO algorithm depends on the physical parameters 
of the IGM, in particular the temperature at the aver- 
age density, the UVB ionisation rate and the index of the 
temperature-density relation. We adopted for these parame- 
ters the values from the cosmological hydro-simulation used 
to validate the FLO performances in paper I. In turn, these 
values are in agree ment with the most r ecent observational 
measurements (e . g. Ricotti et all 2000l: Schave et al 2000l : 
IScott et al]|200d : iTvtler et alj|2004 iBolton et alll2005h . 

In Fig. [7] we showed that values of the UVB ionisation 
rate Fi2 ^ 2 (corresponding to UVB intensities at the Ly- 
man limit, log Jll ^ —21.15) are rejected at more than 3er 
because they over-estimate the average density of the IGM. 

We investigated how the variation of the other two rel- 
evant parameters, To, 4 and a, influence the reconstruction 
of the 5-field with FLO. The results are shown in Fig. [9] A 
temperature at the average density almost a factor of two 
lower than the reference one, which implied a decrease in 
the smoothing scale of ~ 30 percent, could reconcile values 
of Ti2 = 2 with the IGM average density but not values as 
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large as Fi2 = 3. However, this temperature is at the lowest 
end of the measured temperatures in the IGM at the consid- 
ered redshifts. Higher temperatures and high values of Fi2 
allow to recover the correct IGM average density if a lower 
a index is adopted. In particular, an inverted temperature- 
densit y relation with a ~ —0.5, as measured bv lBolton et al.l 
(2008) from the flux probability distribution function, would 
require a high temperature, and a high UVB ionisation rate 
to recover the average IGM density. 

However, the validity of the FLO reconstruction with a 
different set of parameters than the reference one needs to 
be verified with numerical simulations and we plan to do it 
in the next paper of the series. 



5 COMPARISON WITH SIMULATIONS 

We compared the reconstruction of the QSO over-density 
by FLO with the average density field from a cosmological 
hydro-simulation in the proximity of the halos that would 
likely harbour the QSOs in our sample. 

5.1 Simulated lines of sight 

The simulations were run with the parallel hydro-dynamical 
(TreeSPH) code GADGET- 2 based on the conservative 
'entropy-formulation' of SPH fSpringcl 2005]). They con- 
sist of a cosmological volume with periodic-boundary con- 
ditions filled with an equal number of dark matter and gas 
particles. Radiative cooling and heating processes were fol- 
lowed for a primordial mix of hydrogen and helium. We 
assumed a mean UVB produced b y QSOs and galaxies 
as given by lHaardt fc Madaul (|l996l ) with helium heating 
rates multiplied by a factor 3.3 in order to better fit ob- 
servat ional constraints o n the temperature evolution of the 
IGM (Sch ave et al.ll2000l ). This background giv es naturally a 
Tuvb ~ 10 -12 at the redshifts of interest here (|Bolton et al.l 
2005). The star formation criterion is a very simple one that 
converts in collision-less stars all the gas particles whose 
temperature falls below 10 5 K and whose over -density is 
larger than 1000. More details can be found in IViel et al.l 
l|2004l ). The cosmological model corresponds to a 'fiducial' 
ACDM Universe (the B2 series of IViel et aill2004 ). 

We used 2 x 400 3 dark matter and gas particles in a 
120 ft -1 comoving Mpc box. Having a large box size is cru- 
cial since the influence zone of the QSOs is usually of the 
order of some Mpc or tens of Mpc (see Table[T]of this study) . 
The gravitational softening was set to 5 h~ l kpc in comov- 
ing units for all particles. We analysed the output at z = 2.2 
and run a Friend-of-Friend (FoF) algorithm to identify the 
most massive collapsed haloes that should host the QSOs. 
We found about 400 (54) haloes whose total mass is larger 
than 2 x 10 12 M Q h' 1 (10 13 M Q h' 1 ). We then pierced the 
simulated box along the 400 lines of sight (LOSs) intersect- 
ing the centre of the haloes with M > 2 x 10 12 M h' 1 
and along random directions. This latter sample constitutes 
our 'control' sample. We explicitly checked that the correla- 
tion function of the haloes with masses larger than 2 x 10 12 
Mq ft -1 is reasonably well fitted by a power-law function 
with ro = 6 Mpc h T 1 and slope —1.8 i n agreement with ob- 
servational results l|Croom et al.|[2005l ). Furthermore, it was 
recently found that QSOs could be typically hosted in haloes 
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Figure 10. Comparison between the FLO average density field 
for the observed QSO sample (solid dots, upper limit case) in the 
reference model with the F07 slopes and the density field from 
simulations (empty triangles with 1 a error bars). The window in 
the upper right corner is a zoom of the fluctuations at the mean 
level. The shaded area marks the 3 a regions, the dotted lines 
mark the 1 and 2 cr and the dashed one the mean density level. 



of mass 3 x 10 12 Mp h 1 regar dless of their luminosity and 
redshift (|da Angela et al.ll2008l ). 

5.2 Comparison of the density fields 

For each one of the 400 simulated lines of sight, the den- 
sity contrast, the temperature, and the peculiar velocity are 
known pixel by pixel. Peculiar velocities are small, typically 
less than 100 kms -1 , and randomly oriented. However, since 
the observed density fields were recovered in redshift space, 
the mock lines of sight were modified by correcting the red- 
shifts of the density field (z id) with the peculiar velocity 
field to obtain the density field in redshift space (z nC w) using 
the formula w pec (^oid) = c (z nC w — «oid)/(l + (z new ~\~ 2old )/2) 
and the periodic-boundary conditions. Then, the mock LOSs 
were shifted (applying the periodic-boundary condition) in 
order to have the peak of the most massive halo (where the 
QSO should reside) in the first pixel and the proper separa- 
tion from this pixel was computed for all the other pixels. 
The total length of the simulated LOSs is ~ 50 proper Mpc. 
The obtained density fields were finally rebinned into steps 
of 4 proper Mpc as in the case of the observed ones. The 
final simulated density field was computed by averaging in 
each bin the mean <5 of 100 samples of 19 LOSs (the average 
number of objects per bin in the observed sample) extracted 
from the total sample of 400 mock LOSs with a bootstrap 
technique. The error bars on the final density field were com- 
puted from the standard deviation of the 100 samples. 

The comparison between the simulated density field and 
the one recovered from observed spectra is plotted in Fig. 1101 
The simulated <5 field is consistent with the mean IGM den- 
sity down to 4 Mpc and up to 48 Mpc (see zoomed window) . 
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Figure 9. Variation of the FLO output <5-field with the physical parameters of the IGM. The solid line represents the fiducial model 
with To y 4 = 1.8, Ti2 = 1 and a = 0.6 with the F07 QSO slopes. The dotted and dashed lines are obtained for Tig = 2 and IT2 = 3, 
respectively, and varying the temperature and/or the index a as marked above the panels. The shaded regions and the dotted horizontal 
lines are the same of Fig. [7] 



The disagreement between the simulated and observed den- 
sity contrast in the first bin and, in particular, the large value 
of 8 obtained from the mock spectra, could be explained by 
the fact that the considered simulations do not include in 
the most massive halos the presence of the AGN that would 
accrete part of the gas in the peak. 



6 CONCLUSIONS 

In this work, we used a sample of ~ 6300 Lya absorption 
lines obtained from high resolution, high signal-to-noise ra- 
tio spectra of 21 QSOs in the redshift range 2.1 £ z cm ;$ 3.3 
and the FLO algorithm for the reconstruction of baryon den- 
sity fields (described in details in Paper I) to investigate the 
baryon distribution as a function of the distance from the 
QSO. 

The effect of the QSO radiation (supposed to be 
isotropic) dominates over the UV background (UVB) in a 
sphere with radius varying from ~ 3 to 15 Mpc for our sam- 
ple, with an average of ~ 8 Mpc. The increased ionisation 
flux modifies the local line number density causing the so- 
called proximity effect, that was used in the past to estimate 
the intensity of the UVB at the Lyman limit. However, the 
local line number density is determined not only by the ratio 
between the UV flux from the QSO and the UVB but also 



by a variation of the density field due to the fact that QSOs 
likely lie on density peaks. 

The recovery is extremely sensitive to the emission red- 
shift and to the UV flux of the QSOs. We spent a significant 
effort to obtain, both from the literature and from the spec- 
tra at our disposal, the best estimates for the systemic red- 
shifts of the QSOs in our sample in order to reduce the un- 
certainties to less than 100—200 km s _1 . For the 16 QSOs for 
which the redshift was determined with Hp, [O ill] (5 QSOs) 
or with the Ol, Mgn UV emission lines, the average veloc- 
ity difference with respect to the de terminations present in 
the literature (e.g. iKim et ah 1 12004 based mainly on the 
high- ionisation UV lines and Lya ) is ~ —500 kms" 1 . For 
the QSO UV flux, we adopted three different models: two 
with a power law with a fixed slope and one for which each 
QSO has its own slope obtained from the comparison with a 
spectral library. The resulting ionisation rates for the three 
models are within a factor of two. 

We improved the performances of FLO with respect 
of paper I by adopting a smoothing scale ~ 30 percent 
smaller than the Jean s length, following the prescription by 
iGnedin fc Hul { 1998). This contrivance allows us to solve 
the problem of over-estimation of the over-densities we had 
in paper I and to correctly recover the mean density in the 
IGM. 
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We applied the FLO algorithm to each QSO spectrum 
and recovered the sample-averaged <5 field as a function of 
the distance from the emitting sources, in bins of 4 proper 
Mpc. The obtained results are described in the following 
sections. 



6.1 QSO over-density from the proximity effect 

i) If only the UVB ionisation rate is adopted in eq[T]of FLO, 
neglecting the QSO radiation, a decrease of the mean density 
field significant at ~ 2 a is observed within 8 proper Mpc of 
the emitting QSO, confirming the presence of a proximity 
effect. 

ii) When also the QSO ionisation rate is taken into account 
an over-density significant at more than 3<r is recovered at 
proper separations of less than 4 Mpc (first bin) . 

iii) The absolute <5-value of the QSO over-density is uncer- 
tain by an overall factor of ~ 3. A factor of ~ 2 variation 
is due to the lower and upper-limit 5-value derived for the 
empty bins during the field reconstructions: corresponding 
to no absorber or to one absorber with the minimum H I 
column density detectable in our sample, respectively. This 
uncertainty can be improved by using more spectra with a 
larger signal-to-noise ratio. Another factor of ~ 1.1 — 1.5 is 
due to the different assumptions for the QSO continuum in 
the region blue-ward of the Lya emission. 

iv) The uncertainty on the systemic redshift of the QSOs 
used to determine the density distribution in their neigh- 
bourhood is 200 kms -1 corresponding to separations 

,$ 800 proper kpc. As a consequence, it should have neg- 
ligible effects on the 5-value of the peak. 

v) We compared the resulting density field with the aver- 
age density field in the proximity of massive haloes (M > 
2 x 10 12 Mq h -1 ) in a cosmological hydro-simulations. The 
two distributions are in reasonably good agreement. 

vi) There is a significant discrepancy between our results 
and the previous determinations of the matter distribu- 
tion around QSOs, obtained with the optical depth sta tis- 
tics (ODS, iRollinde et al]|2005l : [Cuimaraes et alJl200Tt ) at 
redshifts (z cm ) ~ 2.5 and 3.8. The ODS method recovers 
gaseous over-densities extending to scales as large as ~ 15 
Mpc while our overdensity is limited to a region closer than 
4 proper Mpc from the QSO. We would need to increase our 
sample, in particular with QSOs without associated systems 
in order to study if the brightest objects resides in more ex- 
tended peaks. 

6.2 Constraints on the IGM physical parameters 

i) In the hypothesis of a temperature of the gas at the mean 
density To = 1.8 X 10 4 K and an index of the temperature- 
density relation for the IGM a = 0.6, an UVB ionisation rate 
of Tuvb — 10 -12 s _1 gives the correct IGM mean density. 
On the other hand, Tuvb ~ 2 x 10" 12 and 3 x 10" 12 s" 1 
are excluded at more than 2 and 3 a, respectively, because 
the recovered FLO IGM density field overestimates the mean 
density. 

ii) Values of Tijvb > 10~ 12 s" 1 can be reconciled with 
the correct IGM mean density if different combination of To 
and a are adopted. In particular, an inverted temperature- 
density relation with a ~ —0.5 used in the FLO algorithm 



gives the correct IGM mean density for Fuvb — 3x 10~ s~ 
and To ~ 2.3 x 10 4 K. Such large values of the ionisation rate 
could arise as a consequence o f the He H reionisation at z ~ 3 
and of UVB fluctations (e.g. iMeiksin fc White! 120041 ) . The 
performances of FLO with different set of parameters than 
the reference one have however to be tested with numerical 
simulation. We defer the details of this analysis to a future 
paper. 

6.3 Limiting factors and future developments 

i) Most of the high redshift QSOs have redshifts deter- 
mined from UV emission lines which are known to be 
systematically shifted with respect to systemic redshifts. 
This is particularly critical for proximity effect studies both 
along and transverse to the line of sight. An improve- 
ment in this sense is expected from the new intermediate- 
resolu tion, UV to near-I R spectrograph X-Shooter at the 
VLT (|Vernet et alj|200g| ). that will be operative from the 
first trimester of 2009. 

ii) The temperature of the IGM gas at the mean density is 
highly uncertain. Its best determinations date back to 2000, 
the evidence of a jump in its value at z ~ 3 was marginal and 
should be verified with the present larger samples of high- 
resolution, high-signal-to-noise QSO spectra. Furthermore, 
those temperature estimates are used in many cosmological 
hydro- simulations to re-normalise the ionising background 
intensity. 

iii) The intensity and nature of the UV background 
is another unsolved riddle which should require a new 
observationally-based determination (with the caveat in i) 
since the present simulations are only now starting to have 
the resolution and the physics (e.g. the radiative transfer) 
needed to derive it self-consistently. 
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